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Abstract 


In the framework of renormalon consideration, a role of the anomalous de¬ 
pendence of gluon propagator on the scale p is shown for a large number of 
quark flavours as well as in the nonabelian theory, where the parameter of the 
’’running” coupling constant evolution has a hidden dependence on the scale. 
A phenomenological matching between the modified renormalon and the static 


quark potential is described to fix the renormalon scales. 

Introduction 

In the QCD perturbation theory, an account for the next-to-leading order term over 
is connected with the uncertainty, caused by a choice of the energy scale, determining 
the Q;<j(/i^) value, since the transition to the /i scale leads to the substitution 



where /^o = (lliV — 2nj)/3, N is the number of colours, Uf is the number of quark 
flavours. So, the physical quantity, represented in the form 


r = roas{f/){l + riOs + 0{al)} 
with the given order of accuracy, will get the form 



r 


roQ;s(h^){l + nas -b 0{al)] . 


Thus, the value of ri coefficient for the ag correction depends on the scale of 
determination. In ref. [Q, one offered to fix the /i scale so that the ri coefficient does 
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not contain terms, proportional to (3q. The latter procedure can be understood as one, 
leading to ri to be independent of the number of quark flavours n/ ~ (3q. As was 
shown in refs.|^, such choice of scale has the strict sense in the framework of the 
l/nj expansion, where the a* correction to the gluon propagator is determined by the 
fermion loop contribution into the vacuum polarization. So, the expression for such 
contribution depends on the regularization scheme, and in the next-to-leading order 
one has 


as{fx'^)D{k^,n) 





Gtt 



( 1 ) 


where determines the transversal part of the gluon propagator = 

S°‘^D{k‘^, + kyk\/k‘^). The constant value C is dehned by the renormalization 

scheme, so = —5/3 and Cy = 0 in the so-called 1^-scheme [|I|. 

The summing of {nfag)"' contributions into the gluon propagator leads to the ex¬ 
pression 

/ 9n 9 N as(e^k^) , , 

as{g^)D{k^,g) = -, ( 2 ) 


i.e., in fact, it results in the account for the ’’running” as value in respect to the gluon 
virtuality. Note, the 1/n/ consideration is exact for the abelian theory, where agD is 
the renormalization group invariant. Moreover, the value does not depend 

on the renormalization scheme, and therefore it results in the scheme-independence for 
the expression 

« =_ — _ 


As one can see in the performed consideration, the transition to the nonabelian 
theory is done by the 2n//3 —>• —/do substitution, that is called as the procedure of 
’’naive nonabelization”, giving the correct ’’running” of the coupling constant. 

Further, consider the physical quantity r, for which the hrst order a^-contribution 
is calculated as an integral over the gluon virtuality with the weight F{k‘^) 


^ = J%- • 

Then the account for the {rifasY corrections to the gluon propagator leads to the 
substitution of ’’running” value as{e^k‘^) for as{fi‘^) 





dvr 




( 3 ) 


so that in the offered procedure of the scale hxing |Q, the introduction of the hrst 
order n/Os-correction results in the substitution /i^ —= p^g^pjC -|- {\n k"^ / , 

where the average value is determined by the integral with the F weight. However, in 
some cases P, P| the integration of the n-th order term in expansion (P) leads to the 
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n! factorial growth of the coefficients for the expansion over a”(/i^), so that this rising 
can be determined by the region of low virtualities as well as large ones of gluon. Then 
one says about the infrared and ultraviolet renormalons, respectively ^ . As was 
shown p, P], the renormalon leads to power-like uncertainties of the r evaluation, so 


Ar 



( 4 ) 


where k is positive for the infrared renormalon and it is negative for the ultraviolet 
one. Indeed, consider a quantity JF, presented in the form 


OO 


na.) = Y. 

n=0 



( 5 ) 


and its Borel transformation, defined as 


OO £ 

B\T](n) = hS(u) + 

n=0 


SO that JF can be restored due to the integral operation 

^(a.)=/”*BM(«)exp(-^^), (6) 

In the case of the factorial divergence in the coefficients of expansion (P), the Borel 
transformation has singular points over u P, 0, so that after transform (p), a pole 
contribution at a point Uk = k/2 




ABm{u) = 

U-Uk 

depends on the rule of a way out the singularity, and 

dvTMfc ^ / ,, \ -2u 


AjF ~ ttk exp 






For the two-point correlator of heavy quark vector currents, for instance, one has 
k = 4 and the corresponding uncertainty can be, in fact, eliminated in the procedure 
of definition for the nonperturbative gluon condensate, having the same power over 
the infrared parameter [p|. However, for the renormalized mass of heavy quark the 
infrared renormalon with k = 1 does not correspond to some definite quark-gluon 
condensate, and the value Am(/i) ~ Aqcd can not be adopted into a definition of a 
physical condensate. 

In the present paper we modify the renormalon in the framework of QCD one-loop 
renormalization group and in the l/uf consideration by the account for the anomalous 
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dependence of the gluon propagator on the /i scale, so that under the ’’naive non- 
abelization”, the modihcation appears in the determination of the evolution parameter 
for the ’’running” coupling constant in different schemes of the renormalization. The 
standard renormalon corresponds to the specihc choice of the hxed point of the gluon 
propagator dependence on the scale. The nonabelian evolution of the gluon propagator 
leads to a hidden scale dependence in the renormalon, because of the anomalous de¬ 
pendence on the jj, scale. A phenomenological matching of the renormalization group 
invariant parameters of the modihed renormalon Aqcd and jig, corresponding to the 
singularity points in the QCD coupling constant and in the gluon propagator, respec¬ 
tively, with the tension of string, giving the linear rise of the potential in the heavy 
quarkonium, and with the constant of the subleading coulomb interaction allows one to 
£x the physical parameters, which, in such way, point to a deviation from the standard 
renormalon on the la confidence level. At asymptotically high virtualities of gluon, 
the presence of the second scale in addition to the evolution parameter of the QCD 
’’running” coupling constant is not essential. However, near the infrared region, the 
presence of the additional scale has to be included into the consideration, as it has for 
the potential in the heavy quarkonium or, for instance, for a construction of models 
for such nonperturbative quantities as the gluon condensate. 

In Section 1 we use the one-loop renormalization group of QCD to find the anoma¬ 
lous scale dependence of the gluon propagator. In Section 2 we modify the renormalon 
in the framework of 1/n/ consideration and under the procedure of ’’naive nonabeliza- 
tion” for the ’’running” constant. In Section 3 we find the hidden scale dependence in 
the renormalon, because of the nonabelian evolution of the gluon propagator. In Sec¬ 
tion 4 we make the phenomenological matching of the modihed renormalon parameters 
with the potential in the heavy quarkonium. In Section 5 we use a toy ansatz, moti¬ 
vated by the modihed renormalon, to make a model estimate of the gluon condensate 
value. In Conclusion the obtained results are summarized. 


1 Anomalous dependence of gluon propagator on 
scale 


In the covariant gauge, the gluon propagator has the form 

^ab 


DZik\A) = 




-9v\ + (1 )o;(/i ,k )) 


2 


where aiuj is the renormalization group invariant, since the renormalization group trans¬ 
formations make ’’bare” quantities, independent of the scale, in the form 


(Xi Z^Qji , 

= Z^^ui , 


00 


B 
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where Z3 is the renormalization constant of the gluon held In the 

MS-scheme, one-loop calculations give 

= 1 + i ^(A'(13 - 3a,) - iUf) . 


can be obtained from by the substitution 1/e —> l/e-|-ln47r — 7^;, where 
7 e = 0.5772... The differential equation for ai has the following MS-scheme form 


dai{T) 


a.. 


a)(r)[(1377 - 4n/) - 377a)(r)] , 


dr 12 tt 

where r = Kqcd)) and as(r) = 271/{(3qt). Eq.(|^) has the solution 

(13A^ — Anf)CaT'^°' 


az(^) = 


where 


Ua = 


1 -|- 3iVCar”“ 

13A^ — 4n/ 
6^0 


and Ca is the renormalization group invariant. 

Then eq.(D allows one to represent the gluon propagator in the form 




i:ab 






^{-9uX + {l-b)A^^ , 


( 7 ) 

( 8 ) 


( 9 ) 


where b is the arbitrary gauge parameter^, being the renormalization group invariant. 
In eq.(|D we take into account the anomalous dependence for the gluon propagator on 
the scale, to = only. 

The choice of hxed point a; = 1 corresponds to the standard renormalon. The 
consideration of the gluon propagator at a; 7^ 1 leads to the modihed renormalon. 


2 Modified renormalon in l/uf consideration 


In the leading order over 1/n/ one has Ua = 1, and the cases with 1) a;(/io) > 1, 2) 
a;(/io) = 1 and 3) a;(/ro) < 1 at some large /io values can be formally come to the 
different choices of the renormalization group invariant fig, such that uj{fi’^g) = 0 and 
1) fig < Aqcd, 2) fig = Aqcd and 3) fig > Aqcd, respectively. The a^ifi^) value is 
considered as the formal one-loop expression 




dvr 


Poln 


^^QCD 


^In ref.p the solution with the fixed 6 = 3 value is considered and uj is expressed through 6/a;. 
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( 10 ) 


Thus, the l/uy consideration gives 

a;(/x2) Q(^(/i2) _ 


where d is the scheme-independent invariant of the one-loop renormalization group, 
and it is equal to 

, _ 1 ^QCD _ 47r 

Further, one can repeat the calculation of {rifas)'^ contributions into the gluon propaga¬ 
tor with the substitution of value (^) for so that this replacement corresponds 

to the account for the anomalous dependence of the gluon propagator versus the /i 
scale. Then the scale hxing due to the next-to-leading order correction results in 
the expression 

= /i^ exp{C -I- (In k'^/ /i^) -|- d} , 

and the summing of the corresponding contributions modihes eq.(@) 




as{e^~^'^k‘^) 


asje^k^) 

uj{e^k?‘)k?‘ 


( 11 ) 


Next, one can dehne the Id-scheme, introducing Cy = —d. Then 


aV _ 5/6-d/2AMS 

^^QCD ~ ^ ^^QCD 


In the fd-scheme the perturbative potential between the heavy quark and antiquark in 
the colour-singlet state will have the form 


V(q^) 


4 47ra^(g^) 
3 


( 12 ) 


at = k^. Potential (|^ comes to the Richardson’s potential |^, when one uses 
(g^-|-A^) instead of (g^), and A hxes the linearly rising part of potential, conhning 
quarks with the distance increase, 

AVii„(x) = a|x| = ^ AL|x| . (13) 

The fitting of mass spectra for the charmonium and bottomonium in the Richardson’s 
potential gives 

Alin = 398 MeV . 

However, one usually takes a = 0.18 ± 0.01 GeV^ as the scale normalization for the 
lattice computations, so that this string tension is equal to that of the Cornell model 
[^]. Therefore, we will use 

Alin = 425 ± 15 MeV . (14) 
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Value ( p!^ has to be compared with the evolution scale A of the ’’running” coupling 
constant of QCD, where = 0.117 ± 0.005 corresponds to the one-loop value 

^QCB — 85 ± 25 MeV |§]. If one assumes, that A does not depend on the number 
of flavours, then at the scale of ’’switching on (off)” the additional flavour of quarks 
fi = rrinf+i, the QCD coupling has a discontinuity, related with the step-like change 
l3o{nf) —>• Po{nf + l) = /3o(n/)—2/3. One can avoid such discontinuities, if one supposes 
that A depends on the flavour number, so that as(/i^ = = as(/i^ = 

-f 1). Then in the one-loop approximation one hnds 


Setting A*^^) 


_ /^(nf+1) ( ^nf+1 \ 3/3o(-/) 

VAO/+i)y 

(15) 

85 ± 25 MeV, one gets 


e-d/2 ^ ^ 42 ^ Q_4Q _ 

(16) 


Hence, the rescaling, accounting for eq. ([T5|) , results in the d value, that agrees with 
the standard renormalon {d = 0) within the current accuracy up to Icr, corresponding 
to the conhdence probability, equal to 30%. 

However, the quantity asD, obtained in the framework of l/uf consideration under 
the procedure of substitution 2nf/3 —>• —Po^Uf), possesses the renormalization group 
properties, which disagree with the scale dependence of OgD in the nonabelian theory, 
where Ua % 1. The restoring of correct properties in the nonabelian renormalization 
group for the gluon propagator results in a hidden scale dependence of the renormalon. 


3 Nonabelian generalization and hidden scale de¬ 
pendence 

One has % 1 in the nonabelian theory. Introduce the effective coupling constant 

- / 2 dvr 

where the evolution parameter A(/i) has the following dependence on the /i scale 


A(/i) — Aqcd 


/i 




(17) 


^AqcD' 

Thus, expansions in the perturbative theory of QCD built over the powers of d<j, that 
has the following transformation properties, caused by the hidden scale dependence. 




a,(/i?,A(/i)) 


1 + In ^ 


2 ? 


(18) 
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a^(/i2, A(/X2)) = 




1 + , A(/ii)) In ^ 


A'(/i 2 ,/^l) , 


(19) 


where 


A’(h2,/il) = 




One finds F = 1 at a; = 1 or rza = 1, so that the choice of fixed point or the abelian 
case make us to return to the consideration, performed in the previous section. Then 
A(/i) = /ig, i.e. we come to the Id-scheme, and the standard consideration with = 
Aqcd becomes valid for the choice of fixed point. 

Further, it is evident, that the introduction of the fermion loop contribution into 
the gluon propagator results in the expression 




as{e^, A{fi)) 


k^ 


k^ 


( 20 ) 


where the hidden dependence on the fi scale in A(/r) or in d{fi) restores the correct 
renormalization group properties of the agD quantity in the nonabelian theory, so that 


Accounting for eq. (|^) , one can easily get the generalization of equations, giving the 
BLM-procedure of scale fixing in the ’’running” constant and the renormalon modifi¬ 
cation under d —>• d(/r). 

However, physically measurable quantities have to be invariant under the renormal¬ 
ization group transformations, so that corresponding expressions can contain the agD 
value, given at a physical scale, being invariant under the renormalization group and 
defining a normalization condition. 


4 Phenomenology of scale fixing 

Let us consider the momentum-space potential in the heavy quarkonium at 
the Id-scheme under the Richardson’s prescription 

v( = — 

’ 3 g2 o;(/z2(g2)) ’ 

where /i^(g^) = + fAi, so that the expansion over the /fAi powers gives 


0 in 


V{q^) = 


4 dvr /Svr 47r Oisi/jA) 




3 g2 \pQ g2 pQ 


+ 0 ( 1 ) , 


2 


( 21 ) 










so that at large |x| values and at /3o(n/ = 3) = 9 one gets 


y(x) 


Idvr 2 i I 4 1 /dvr 

IfwlT 



( 22 ) 


Note, that the coefficient at the scale, giving the cr|x| term in the potential, is twice 
greater, than in eq.(^). If one accepts the given scheme of matching of the gluon 
propagator parameters with the parameters of potential, conhning the quarks in the 
heavy quarkonium, then at /i^ = /ig(l + 0{q ‘^/one will have 


A(/r) ~ A{fig) = fig , 


and we obtain, that the effective ’’running” coupling constant of QCD is given in the 
I/-scheme. Moreover, the Ann value, determined in Section 2 in the l/nj consideration, 
is related with the evolution parameter of the coupling constant by the expression 


™ = v/2 , 
Ag 


that is in a good agreement with the experimental value in eq.([l^). Thus, the matching 
of the evolution scale fig with the gluon string tension in the heavy quarkonium results 

in the value _ 

fi^^ = 136 ± 4 MeV (23) 

at Uf = 3. 

Further, as was shown, the account for the fermion loop contribution into the gluon 
propagator gives the potential of the Richardson’s form, with A(/i) = fig- In such 
potential, the expansion over the /fi^ powers gives 




4 dvr /dvr fig 27i\ 
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+ ^ + 0 ( 1 ). 


(24) 


Comparing eg .(^11) and eq.(^4|), one notes, that the modihed potential (^) at large 
distances can be obtained from the Richardson’s potential by the following 


V'(x) « ^ + V 


Rich ( 


(25) 


if one supposes 


i.e. 


( 2n _ 47r 

^s\Ag) ^ 5 


In 


fi 


9 


^QCD 


= 1 . 


(26) 
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V^(|x|), GeV 


1 

0 

-1 

-2 

fm 

Figure 1: The potentials, corresponding to the Cornell model (dots), the model with 
the ’’running” asjio value (solid circles) and with the ’’running” as value plus the linear 
correction (empty circles) 



The additional term in eq.(^) appears as the correction to the linear potential, conhn- 
ing the quarks, and its presence is due to the modihcation of the renormalon. Further¬ 
more, the matching of the modihed potential with the Richardson’s potential over the 
coulomb deviation from the linear rise at large |x| allows one to £x the ratio ^g /^qcd 
or asi^"^) (see eq.(^^). As has to be expected, the presence of the additional scale 
in the agD quantity becomes unessential at low distances, i.e. at large momentum 
transfers, when the dependence on the scale is determined by the coupling constant 
evolution. 

The comparison of the modihed potential with eq. (^Sf) and with the Cornell model 


potential is given on hgure 1. As one can see, at large |x| > 0.5 fm, where the 
modihcation is valid, the Richardson’s potential with /r^, determined in eq. (p^) , and 
with the account for the additional linear term is the quite accurate approximation, that 
is also close to the potential of Cornell model. As was noted 0, the QCD-motivated 
potentials, fitting the spectra of charmonium and bottomonium, give one and the same 
form of the R(x) dependence at 0.2 < |x| < 1 fm with the accuracy up to an additive 
shift AV, independent of |x|. This form is caused by the transition of the coulomb¬ 
like potential, following from the asymptotic freedom of QCD at low distances, to the 
linearly rising string-like potential, conhning the quarks at large distances |x| > 1 fm. 
In the present section we have found the connection between the scale of the QCD 
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running” constant evolution and the string tension. 


5 Model-dependent estimate of scales 


A toy ansatz, based on the extraction of the pole singularity of u in the gluon propa¬ 
gator in the euclidian space and in the Id-scheme is given by the expression 


Dmodik^) 


iln^fig/AQCp) 


that is used for the one-loop calculation of the gluon condensate {—G^) value. To get 
the renormalization group invariant expression in the dimensional regularization, one 
has to extend Dmod into the (4 —2e)-dimensional space, so that = Dmod{^ + ^Zi): 

where zi is a single-fold determined value, because of the requirement of the renormal¬ 
ization group invariance of the gluon condensate value. Straightforward calculations 
give 



vr2 Aqcd 



^QCD 




(27) 


Note, that in accordance with eq.(|2^), the gluon condensate tends to zero, if one has the 
purely perturbative-like expression for D^od at = 0 or if one takes the hxed point 
fig = Aqcd- The substitution of numbers for the scales, estimated in the previous 
section, results in 

{—G^) ^ 0.011 GeV^ , 

TT 


(28) 


where fig = e^'^fig = 316 ± 12 MeV. Value (^) is in a good agreement with the gluon 
condensate estimate in the framework of QCD sum rules . 


Conclusion 

In the present paper we have shown, that in the l/uj consideration, the account for the 
anomalous scale dependence of the gluon propagator leads to the modihcation of the 
renormalon, because of the redehnition of the evolution scale for the ’’running” coupling 
constants of QCD. When one restores the correct renormalization group dependence of 
the gluon propagator on the scale in the nonabelian theory, the modihed renormalon 
obtains the hidden scale dependence, that is expressed in the specihc transformation 
properties of the effective ’’running” coupling constant. The phenomenological match¬ 
ing of the modihed renormalon with the tension of gluon string between the static 
quarks and with the subleading coulomb interaction at large distances results in the 
hxing of the two renormalization group invariant scales of the modihed renormalon, 
and it gives the normalization condition for the ’’running” coupling constant. The 
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modification of renormalon leads to the linear term, additional to the Richardson’s 
potential. Furthermore, we have made the model estimate of the gluon condensate 
value at the presence of two infrared parameters, determining the scales of the pole 
singularities in the ’’running” QCD constant and in the gluon propagator up to one 
loop. At large virtualities of the gluon, the coupling constant is given by the expression 




27r 

/?oln(m//ig) ’ 


to one loop, so that 


= 3) = 136±4MeV, 


that is in agreement with the experimental estimates of Q!s(m|). Moreover, 

271 


(/^g) 


/3o\n{fig/A qcd) 
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